Spreading of Block Copolymer Films and Domain Alignment at Moving Terrace 
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We investigate spreading of phase separated copolymer films, where domain walls and thickness 
steps influence polymer flow. We show that at early stages of spreading its rate is determined by 
slow activated flow at terrace steps (i.e. thickness steps). At late stages of spreading, on the other 
hand, the rate is determined by the flow along terraces, with diffusion-like time dependence i -1 / 2 . 
This dependence is similar to de Gennes and Cazabat's prediction for generic layered liquids-, as 
opposed to the classical Tanner's law of drop spreading. We also argue that chain hopping at the 
spreading terrace steps should lead to the formation of aligned, defect-free domain patterns on the 
growing terraces. 

PACS numbers: 81.16.Rf, 82.35.Jk, 47.54.+r, 83.80.Uv, 68.15.+e 



I. INTRODUCTION 

Equilibrium structure of block copolymer systems 
was in the center of theoretical and experimental re- 
search for a number of years now, and is mostly 
understood2*^^iSiL2iS. It is known that in the bulk, 
at temperatures below order-disorder transition (ODT), 
molten diblock copolymers phase-separate into periodic 
structures with periodicity dependent on polymer size. 
The actual structures formed depend upon the propor- 
tion of different blocks inside copolymer, and include 
lamellar, cylindrical, spherical and various bicontinuous 
morphologies. 

The spreading dynamics of copolymer Alms, on the 
other hand, is much less investigated^. While little dif- 
ference besides significantly higher viscosity is expected 
in the spreading of macroscopic copolymer drop, large 
deviations may arise on microscopic length scales, when 
the film thickness approaches domain period. In particu- 
lar, thin films inherit domain structure of the bulk sam- 
ples, but their domains are typically aligned parallel or 
perpendicular to the substrate. Combined with extreme 
anisotropy of diffusion relative to the domain orientation, 
spreading dynamics of thin films is expected to be very 
different from that of simple fluids. 

In the present paper we look at the spreading of phase- 
separated copolymer films on wetting substrates in the 
regime of complete wetting. It is known that copoly- 
mer films at small thicknesses adopt a terraced topology, 
with thickness of each terrace defined by the underlying 
domain structure 12 . The copolymer domains are then 
localized inside these terraces (Fig. [TJ, so that terraces 
can slide on top of each other. 



*The following article has been submitted to the Journal of 
Chemical Physics. After it is published, it will be found at 
|http://ojps.aip.org/jcpo| 



General spreading of layered drops was first investi- 
gated by de Gennes and Cazaba1ii*i£, and consists of vis- 
cous flow along each terrace, combined with disjoining 
pressure driven flow at terrace steps. Here we investi- 
gate spreading dynamics of the layered copolymer films. 
On one hand, the thickness of copolymer terraces (layers) 
significantly exceeds atomic sizes, so that effect of disjoin- 
ing pressure is negligible at all but the bottommost ter- 
race. Hence, the flow between terraces is driven by other 
factors, presumably by jump in the internal stress. Ad- 
ditionally, the thickness of each copolymer terrace needs 
not be uniform. This results in the forced diffusive flow of 
material along the terrace with diffusion constant deter- 
mined by the underlying domain topology and strength 
of copolymer block interaction. We find that at large 
times the spreading rate v for the copolymer film takes 
the diffusion like time dependence v cx i -1 / 2 , as in gen- 
eral layered drops of Refs. Q and El and contrary to the 
Tanner's law for simple liquids, which predicts v oc t~ 9 . 

The activated flow of chains at terrace steps, with 
chains hopping between copolymer domains, can also 
account for the growth of aligned domain patterns in 
spreading films. Specifically, asymmetric copolymers 
that phase separate into cylindrical domains in bulk, of- 
ten form striped pattern in the thin films, similar to those 
shown on Figure ^uLiii. Easy control over periodicity of 
these patterns has driven large interest in their forma- 
tion, but lack of ordering and generally large number of 
defects limited their applicabilitj*Si2iS. 

However, Hahm and Sibener— recently reported that 
nearly perfect alignment is observed on annuli struc- 
tures of polystyrene-block-polymethylmethacrylate (PS- 
PMMA). In that experiment, a thin film of PS-PMMA 
was spin cast onto a Si02 substrate that was initially 
covered with a minor polar solvent. Dewetting proper- 
ties of the solvent resulted in the formation of bare circu- 
lar regions with thick annular rims (Fig. EJ). On subse- 
quent annealing the thick rims spread out and flatten. 
The expected striped domains showed unexpected ra- 
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FIG. 1: Sketch of the layered structure of the cylinder forming 
PS-PMMA copolymer discussed in the text. Shaded regions 
correspond to PMMA domains. 



dial alignment. This strong alignment along the spread- 
ing direction was confirmed in later experiments with 
a different geometry— Below we will argue that hop- 
ping of copolymer chains between domains at terrace 
steps, necessary for the spreading, may be responsible 
for this alignment. In particular, chain hopping should 
lead to the continuous growth of cylindrical domains, 
while thickness gradient near terrace steps should gen- 
erate alignment of the growing domains. As a result, one 
may expect aligned, defect-free domain patterns along 
the whole area of spreading, similar to the observations 
by Hahm and Sibeneii 

Improving the ordering of copolymer film domains is 
regarded as a key to many potential application o 15 i 16 i 17 . 
Accordingly, much attention has been given to poten- 
tial alignment mechanisms: local alignment near a ter- 
race stepi^, thermal annealing^, electric fields^, growth 
from an aligned boundarjiiAiiS and imposed shear—. 
The alignment under spreading observed by Hahm and 
Sibener represents a potentially new mechanism, differ- 
ent from these others. To understand this mechanism is 
important for the overall goal of controlling copolymer 
film patterns. 

This paper is organized as follows. In section [H] we 
review the known facts about homopolymer spreading. 
In sections IlIII - IVl we develop our model for the copoly- 
mer film spreading. Then, in section IVII our model is 
compared with the experimental results on PS-PMMA, 
and a new dynamic mechanism of domain alignment is 
discussed. 



II. SPREADING OF HOMOPOLYMER FILMS 

We will start with a brief review of the spread- 
ing dynamics of regular homopolymers. This spread- 
ing has been widely investigated and is mostly 
understood2L22i2ii2i. The rate of this spreading, as well 
as the shape of the advancing liquid front, can be derived 
from the lubrication approximation. Within this approx- 
imation, the two-dimensional current density of material 




FIG. 2: Sketch of annular structures observed by Hahm and 
Sibener, based on Figs. 2 and 3 of Ref. Nearly perfect 
alignment of cylindrical domains in the radial direction was 
observed on the top surfaces of the terraces. In the actual 
experiment terraces extension in radial direction is 50 times 
larger than their thickness. 

flow J is given by 

J =- h i{ b+ l h )^ w 

where r\ is the viscosity, h is drop thickness, b is the 
slip length and p is pressure inside the drop. There 
are generally two major contributions to this pressure: 
one is from Young-Laplace capillary effect 7V 2 h and the 
other is from Van-der-Waals disjoining pressure W(h) = 
—A/(6irh 3 ). The coefficient A is a material dependent 
constant. Thus 

p= -jV 2 h + W(h). (2) 

In the frame moving with velocity U of the contact 
line, the current becomes J' = J — Uh. For steady-state 
spreading, this current must vanish, and the spreading 
profile h{x) is determined by 

V [- 7 V 2 /i + T^(/i)] ^(b+^hj =-U. (3) 

Based on the thickness h of the drop, there are three 
distinct regions that can be distinguished. First of all, 
there is a region of macroscopic drop where film thickness 
is much larger than the slip length (ft> 6), and contri- 
bution from disjoining pressure is negligible. Equation 
P)l in this region reduces to the well known spreading 
equation investigated by Tanner— 

^-h 2 V 3 h = U (4) 

which at large h predicts zero curvature wedge-like profile 
h cx |x| In 1 / 3 \x\ and cubic dependence of the spreading 
rate on the apparent contact angle U cx -^O 3 . 

The second distinct region of is called the precursor 
film*^, i.e. the area of the film where thickness is very 



3 



small so that plug flow takes place and Van-der-Waals 
interactions dominate. Then @ reduces to 



^VW(h) = ^-h~ 3 Vh = -u, 
r\ ztttj 



(5) 



which predicts h oc |a;| _1 ' 2 thickness profile. 

Finally, there's an intermediate region characteristic 
only to polymeric drops and denoted as the foot region. 
Existence of this region is associated with very large slip 
length of polymeric liquids, so that there is a region of 
the film where Van-der-Waals interactions are already 
insignificant, but flow is still of the "plug" type. Then 
(PJ reduces to 
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KS7 6 h = U, 



(6) 
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with film thickness obeying 2 ^ h oc 

The described model for the polymer spreading suf- 
ficiently well describes both precursor and macroscopic 
film regions. However, a; 3 / 2 dependence for the foot is 
not generally observed^ 3 .. 

Similar effects may be expected in copolymer drops. In 
the precursor film, for example, Van-der-Waals interac- 
tions with the substrates significantly exceed repulsive in- 
teractions between copolymer blocks, so that copolymer 
should obey the same dynamics JSJ as a homopolymer. 
Similarly, in the macroscopic drop region, domain orien- 
tation may be expected to be random and macroscopic 
viscosity remains isotropic. Hence, the domain structure 
may be accounted for via a macroscopic viscosity so that 
general spreading dynamics still follows the Tanner's law 
©■ 

However, in the foot region, one may expect even 
stronger deviations from than those in homopoly- 
meric liquids, as periodic domain structure of copolymer 
gets aligned parallel to the substrate and a quantized 
thickness profile emerges. Then the continuous descrip- 
tion of equation becomes invalid. To describe the 
flow of material inside this terraced structure we will 
look at two different regions of material flow: on one 
hand, there's continuous flow inside each terrace, which 
can be described from a hydrodynamical standpoint; and, 
on the other hand, there's activated flow at the terrace 
steps, where chains have to hop between discontinuous 
domains of different terraces. 



III. FLOW INSIDE A SINGLE TERRACE 

Hydrodynamic Approximation. The terraced structure 
of copolymer films comes as a result of domain alignment 
parallel or perpendicular to the substrate. Some of these 
terraces may slide on top of each other (Fig. [IJ, with 
flow of material inside each terrace being strictly one or 
two dimensional and parallel to the substrate. Then it is 



convenient to characterize profile of each terrace in terms 
of area density p oc h, so that the flow can be described 
by hydrodynamic equations of motion: 



dp 
dt 



V ■ (pu) = 0, 



(J) 



and 



Dt 



Vi 0~i 



kllj 



(8) 



where k is the viscous friction opposing the flow, and cr^ 
is the stress tensor. In the discussion that follows, we 
will consider only slow (Du/Dt sa 0) spreading in the 
x-direction. Additionally, we will neglect y-dependence 
of the stress tensor a. This dependence, which arises 
from the domain-induced anisotropy of phase-separated 
copolymer and potentially affects such parameters as 
direction of spreading, is insignificant for the present 
discussion and will be briefly considered in section IVII 
Hence, (jSJ reduces to the familiar Darcy's law 
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da xx da x 



dx 



dy 



1 da xx 
k dx 



(9) 



Non-uniform stress o~ xx exists in the film because it is 
not at its equilibrium thickness fiQ. The stress can be 
expressed in terms of the free energy per chain F c . This 
free energy evidently has a minimum when h = ho, and 
for h near ho F c is quadratic in (h — ho) 
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F c m F + go^fABVo 



(h-hpf 

h 3 
"o 



(10) 



Here vo is the volume displaced by each chain. We 
treat the polymers as an incompressible liquid; thus vq is 
a fixed constant. The jab is the interfacial energy of the 
interface between the A and B blocks. The quantity g is 
a numerical constant that depends on the domain mor- 
phology and how it is distorted in the process of changing 
h (e.g. go — 3/2 for plane lamellar terraces 2 ^). 

To determine the stresses, we consider a small section 
of the terrace of width Ax and length Ay with thick- 
ness h ^ ho- The section is prevented from relaxing to 
height ho by forces from the adjacent film. If the section 
were detached from adjacent terrace, it would contract 
or expand laterally, allowing the height to return to its 
equilibrium value. To prevent this, the adjacent terrace 
must exert a tensile stress on the section. In particular, 
during spreading in x direction, Ay stays constant and 
stress a xx is 
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dF T 



hAy d{Ax) 



4 



Using incompressibility of polymer hAxAy = V and 
the fact that the terrace's free energy Ft = (V/vq)F c we 
finally get 



h dF c 
Vq dh 



9 jh-hp) 

~2g To — 1AB 



h 2 
"o 



(11) 



The numerical constant g depends on the go of Eq. 1101 
and in the trivial case of lamellar domain g — go- I n 
general, the tensile stress also includes a Laplace- Young 
contribution j'S/ 2 h from the curvature of the free sur- 
face. This part of the pressure is negligible in the long, 
flat terraces we are considering. Then the profile of the 
spreading terrace is given by 



diffuse between them with velocity Ud- The net velocity 
of material flow is then 



tu s + (1 - r)u d , 



(14) 



where r is the fraction of the sliding chains, i.e. chains 
that are in contact with the substrate. If the fraction r 
of sliding chains is small, we can temporarily switch to 
the frame moving with velocity Ud where sliding chains 
become diffusing, with the drag force on them given by 
Einstein relation /q = £ji^-(u s — u d ). Here D c is the 
self-diffusion coefficient for chains moving along domain 
walls. The net force between sliding and diffusing chains 
is then / = Tfo/v = r^-(u s - u d ), and the Darcy's 
law for each of these components is 



dh 
~di 



d_ 

dx 



ZgjAB h dh 



Kq dx 



(12) 



Away from the terrace steps, the thickness of the ter- 
race is close to /loi so that to the first order in h — ho the 
thickness profile is given by 



(1 - t)V ■ a + Tk d (u d - u a ) =0 (15) 
rV • a + £fc s u s - Tkd(ud - u s ) = (16) 

where kd — k B T '/vqD c , k s — rj/bh and vq is the volume 
of one chain. These equations, together with l|14|) give 
the final expression for the velocity: 



dh 
~di 



D, 



d 2 h 

dx 2 ' 



(13) 



where D h = 2g^ AB /h k h=ho . 

Dissipative Forces at the Substrate. We now consider 
the physical origin of the dissipative coefficient k. The 
most common type of viscous forces are the frictional 
forces at the substrate interface. These forces are most 
easily expressed in terms of the slip length b, i.e. the 
distance below the substrate at which velocity of liquid 
u(z) extrapolates to zero: b = [u/(du/dz)] z=0 — T]/j3. 
Here r\ is liquid viscosity, and (3 is friction coefficient at 
the substrate. The magnitude of the friction force does 
not depend on terrace thickness, so that coefficient k en- 
tering (JHJ is k = k s — (3/ h = r\jbh. 

While this expression for k s is believed to be true for 
macroscopic drops, it fails in thin films due to rapid de- 
crease in monomer mobility near a substrate. For in- 
stance, near a high-energy substrate monomer mobility 
may drop 10 4 times compared to its bulk value 22 , leading 
to a similar drop in k. Introducing phenomenological pa- 
rameter £ = ^buik/ ^substrate, which relates bulk {nbuik) 
and substrate {^substrate) mobilities, the friction coeffi- 
cient for the sliding chains becomes k = £fc s . 

This decrease in mobility falls off rapidly away from 
the substrate. Then the chains that are affected most 
are those in immediate contact with the substrate, while 
remaining chains can freely diffuse between them. This 
leads to a two-liquid model of flow, first introduced for 
polymer spreading by Bruinsma 22 . According to this 
model, the first liquid consists of chains that come in 
direct contact with the substrate and slowly slide with 
velocity u s , while the other liquid consists of chains that 
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(17) 



As expected, this expression reduces to (0 in the case 
when sliding dominates diffusive flow (£fc s <C kd). With 
velocity l|17|) . one-dimensional film spreading is still de- 
scribed by (jT3|l . with collective diffusion coefficient Dh 
being 



Dslid = 



37 



AB 



(1-r) 2 1 



1 



(18) 



Dissipative Forces Between Sliding Terraces. Along 
with dissipation at the substrate, an alternative dissi- 
pation mechanism arises when terraces slide over each 
other. During this sliding motion, brushes from adja- 
cent terraces interpenetrate and viscous dissipation takes 
place. The areal density of the dissipated power TE is 
given by 



I) *"4 US) 



where v is spreading velocity, n* is the effective viscosity 
of one of the copolymer blocks, and £ is the characteristic 
depth of interpenetration of two adjacent brushes. In 
the plane lamellar morphology the brush thickness hj\ oc 
nV 3 so that the interpenetration depth is21: 



C« (R\/h A ) 1/3 ^aN 



.4/9 



(20) 



where Ra is radius of gyration of a free chain, and a 
is a microscopic length scale close to the monomer size 
a i 2 ?' 3 ^ a ~ (kT /^ABalY^ao. The effective viscosity 
rf entering (|19|l is qualitatively different from the ho- 
mopolymer melt viscosity. This is due to the fact that 
chains in a segregated copolymer are bound by domain 
walls and may no longer undergo reptation motion char- 
acteristic of linear homopolymers. Instead, the copoly- 
mer system becomes increasingly similar to the star poly- 
mer melt, where chains are also constrained at their 
vertes^l. The characteristic mode of motion of star poly- 
mers is repetitive contractions and extensions of one arm 
of the star, with time scales exponentially increasing with 
the length of this arm. Correspondingly, the viscos- 
ity of star polymers also exhibits exponential increase 2 ^ 
Tf* w r] arm e Narm / N '' , where N e is the average entangle- 
ment length, N arm = 2£Na/Ha ~ 2N 7 / 9 is average arm 
length, and r\ arm is the viscosity of the corresponding lin- 

7 /9 

ear polymer of length 2N / . Thus the friction coefficient 
between terraces becomes 

= ^ « ^N A ^ ex p (2NT/N e ) , (21) 

so that final dynamics of the terrace is described by 
with collective diffusion coefficient Dh given by 

n ZgiAB GaiAB AT <l/9 ( OM 7/9/ M \ 
Attsc = a ~ A ex P ( ~ 2N A / N e ) • 

(22) 
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FIG. 3: Schematic top view of a step between terraces with 
lying-down cylinders. As thiner terrace (A) grows, material 
from the thicker terrace (C) gets squeezed inside transitional 
region (B), resulting in chain hopping from squeezed bulk 
domains (B) towards stretched surface domains (A) . For sim- 
plicity, only the topmost terraces of the film are shown. 

is the self-diffusion coefficient of a homopolymer melt 
with chains assembled from majority block monomers 
and having the same polymerization as the copolymer 
chain. 

There are two distinct contributions to the above ac- 
tivation energy: the energy cost AEq for the hopping 
block to pass through the incompatible stationary block, 
slightly decreased by the excess energy, or chemical po- 
tential, of the chains inside domains. In stationary sys- 
tems, chemical potentials for all domains are the same 
and no net flux occurs, while in the spreading systems 
chemical potentials may differ between terraces, produc- 
ing net flow of chains. 

To estimate this flow, we notice that the average rate 
at which chains hop between domains with chemical po- 
tential fi and spacing ho is 32 



IV. FLOW AT THE TERRACE STEPS 

Besides continuous flow of material inside each of the 
terraces, flow of material between terraces is needed 
to sustain spreading. However, topologically differ- 
ent domains of adjacent terraces (Fig. are ei- 
ther disconnected, or connected via some intermediate 
morphology* 2 ** In either case, activated events are needed 
to create flow at the step: in the former case, single 
chain hopping is necessary, while, maintaining interme- 
diate morphology in a dynamic system, such as spread- 
ing system of Fig. |3|with step moving to the right, re- 
quires many chain activation events. As the time scales 
for many chains events are significantly longer than for 
single chain hopping, from here on we will assume that 
single chain hopping dominates flow at terrace steps and 
intermediate morphologies do not form. 

Qualitatively, the decrease in chain self-diffusion due 
to activation events has the form^fl* 3 ^ 

n ~ n -AE/kT 

where AE is the activation free energy for the minor- 
ity block to hop through the stationary block, and Dq 



(At)-* = ±£ « ^(-AKa+^T (23) 

For two parallel lamellar layers separated by ho the 
current density jo of chains hopping from one layer to the 
other is given by jo = po(At)~ 1 , where po is equilibrium 
number of chains per unit area. At equilibrium, an equal 
current density flows in the opposite direction. However, 
in our spreading terrace, the thin, leading terrace is not 
in equilibrium with the thicker terrace behind the step. 
The two regions have different values of p, and a net 
current is expected. 

The relationship between the emitting and receiving 
regions may take two general forms, as depicted in Fig. 
4. In this figure the regions of different chemical poten- 
tial are labeled 1 and 2. In Fig. 0Ji the regions lie one 
above the other; in Fig. ^Js, one lies ahead of the step, 
the other behind it. In view of Ea l23l the current of poly- 
mers must be largest in regions with smallest activation 
free energy AE and largest Ap. In the configuration 0Jd, 
the change of /j is concentrated at the step and, hence, 
so is the current. Similarly, in configuration 0Ji domain 
distortion is maximum near the step and the free energy 
of the transition state may be reduced. For example, the 
hopping distance for a B block may be reduced, and thus 
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a) 




FIG. 4: Two cases of stress propagation between terraces: a) 
if domains inside top (2) and bottom (1) terraces are sepa- 
rated, stress from pulled bottom terrace does not propagate 
into top terrace; b) if two terraces are linked (2), stress from 
the pulled terrace (1) equally splits between the two linked 
terraces (2). 
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FIG. 5: Profile of the advancing terrace, 
and net current l|24l) between terraces becomes 

i « -2j%„^^ e (-^^)/ ff , (26) 

where g is between go and 2go depending upon the inter- 
action between terraces. 



the amount of stretching required for a B block to make 
the transition should be reduced^. Thus we expect the 
current to be concentrated at the step in both situations. 
We shall assume that this is the case in what follows. 

The length of the step region where the current is con- 
centrated has a length in the x direction of order ho- The 
current i per unit width of the step is thus approximately 
the net current density times ho'. 



ho A j 



2D p e 



-AEg/kT 



ol i h /kT _ ut/kT 



(24) 



The chemical potential for a terrace of width Ace and 
length Ay to the first order in h — ho is given by 



dF T 



dn 



Mo + 2^0 



h — ho 



-lAB, 



(25) 



Ax, Ay 



where = JAB^o/ho is interfacial energy normalized 
to one chain, Ft — hAxAyF c /vo is terrace free energy 
(see Eqs. fTOl- HT)l. and n — hAxAv/vo- If terraces can 
slide on top of each other (as in Fig. , then the spread- 
ing pressure from the lower terrace does not penetrate 
into the higher terraces, thus leaving thicknesses and 
chemical potentials of the latter at equilibrium. Then 
A M = 1gi* AB (h-ho)/ho- 

On the other hand, if two terrace are bound by domain 
walls (as in Fig. [3] and Fig. QJ)), then spreading pressure 
penetrates past the terrace step. From the balance of 
forces at the step, the thickness h of the terrace on the 
left (Fig. HJd), and thickness h 2 of the double terrace on 
the right are related via (h — ho) = 2(/i2 — ho). Then 
Am = gj* AB {h- h )/h . 

In either case, difference in chemical potentials is pro- 
portional to the deformation h — ho of the lower terrace, 



V. SPREADING OF THE COPOLYMER FILMS 

It was shown above that the spreading of copolymer 
films consists of two parts: flow of material inside each 
terrace, governed by diffusion law (|13fl . and flow of ma- 
terial between terraces, with the flow rate described by 
(|26|l . At early stages of the spreading, when terrace lat- 
eral extension is small, one may expect the dominant con- 
tribution to the total rate to come from the flow at terrace 
steps, while at late times, the dominant role should shift 
to the flow inside terraces. 

To find this time dependence, we turn to the spreading 
of the lowest terrace. The advancing front of this terrace 
(x = L(t), Fig. [HJ is assumed to be pulled with a con- 
stant pressure from the wetting forces. As pressure and 
thickness of a terrace are proportional (eq. JUJ), this 
is equivalent to stating that advancing front is held at a 
constant thickness hf. 



L{t) 



(27) 



The boundary condition at the step between terraces 
(x = 0, Fig. |SJ arises from the balance between the cur- 
rent i|2rj|) entering from the higher terrace, and the diffu- 
sion current i = —^P-D§p- inside the terrace, where D is 

ho ox " 

collective diffusion constant entering Hlri[) and defined by 
(HU or Hence 



-d~x +A{h 



ho) = 0, atx = 0, 



where 



2 -A)7^J_ e (- 
y D kT h 



AE + W )/tT 



(28) 



(29) 
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Finally, the lateral extension L of the terrace is re- 
lated to the current i t at the advancing front via it = 
Po(hf/ho)(dL/dt). Hence 



dL 
~dt 



D_ dh 

hf dx 



(30) 



--L(t) 



with 



L\t=o = 0. 



(31) 

Thus we come to a system of equations ifHH 1271 1251 15171 

I31|l . Exact solution of this system is not trivial, so we 
will look only on large and small time asymptotics. 

The large time behavior is best analyzed in new coor- 
dinate z — x/L(t), which eliminates the moving bound- 
ary condition 1)2 7[l. Differential equation (|13|l with mov- 
ing boundary condition, similarly to the case with sta- 
tionary boundaries, has only one solution^. Then, in 
the z coordinate, the system is likely to asymptoti- 
cally approach some stable state hi(z). Assuming that 
h(z 7 t) = h l (z)-h 2 {z)t-l 3 , L = at a , we get a = (3 = 1/2, 
and 



hi(z) = h f + (h - hf) 
exp (—R 2 z 2 



erfifa 



h 2 (z) = {ho - h f )- 



eiiR 
erfii?z 



crfii? 



(32) 
(33) 
(34) 



hn—hf 



where R — , 
times is given by 



Thus velocity of spreading at late 



>D ho 



V ,-1/2 



(35) 



The early times behavior, on the other hand, is charac- 
terized by nearly constant thickness h w hf of the grow- 
ing terrace (Fig. |SJ. Then all the material entering a 
terrace (eqs. 1261 and l2*5j) should go towards lateral exten- 
sion of that terrace with spreading velocity given by 



dL 



h 



carly dt Pohf 



(36) 



hohf kT 

The characteristic cross-over time between the two ve- 
locities is found by setting w e ariy ~ fiate : 



VI. SPREADING OF PS-PMMA FILMS ON SI0 2 
AND DOMAIN ALIGNMENT 

Spreading rate of PS-PMMA. The theory above pre- 
dicts the absolute spreading rates and thickness profile 
to be expected in the experiments of Ref 0- We may 
thus test our theory by comparing these expected fea- 
tures with the observed ones. In these experiments, a 
thin film of 80,000 amu PS-PMMA was spin cast onto 
a silicon substrate that was initially covered with a mi- 
nor polar solvent. The presence of the solvent resulted in 
the formation of annuli structures, likely caused by the 
dewetting of the substrate, and consisting of a tall rim, 
bare substrate inside the rim, and uniform film outside 
(Fig. [2J. Later, as the solvent evaporated, the material 
started flowing downwards onto the bare substrate area 
inside the rim. This flow is presumed to be driven by 
wetting forces, so that its rate should be expressed by 

The structure of asymmetric, cylinder-forming PS- 
PMMA copolymers on SiC>2 used in this experiment is 
well known^iiS, and consists of plane lamellar layer on 
the substrate, followed by a series of layers (terraces) 
with cylindrical domains inside (Figure To estimate 
the spreading rate of this system, we turn to the bottom- 
most terrace pulled by the wetting forces. The collec- 
tive diffusion constant for material flow inside this ter- 
race is given by (|18l) and contains two separate contri- 
butions. On one hand, there's a friction on the terrace 
as a whole (second term in Eq. I18fl . However, the hy- 
drogen bonding between PMMA and SiC>2 is very strong 
(£ 3> 1), and contribution from this term to the dif- 
fusion constant is small. The first term, on the other 
hand, describes self-diffusion of fast chains through the 
matrix of substrate-bound slow chains. Since there are no 
walls obstructing this flow in the bottommost terrace of 
PS-PMMA, contribution from this term dominates and 
Dsiid ~ 3(1 - T) 2 v Q D c j AB /k B TTh . 

The fraction r of chains that are in contact with the 
substrate can be estimated based on a strongly-stretched 
approximation^, which predicts that the density of chain 
ends in a flat brush is g[r) oc r/(h 2 — r 2 ) 1 / 2 . Since 
strongly stretched approximation neglects deviations of 
chain paths from their mean, and all paths are pre- 
sumed to go unidirectionally towards the free end of the 
brush, we can assume that all chains ending within a 
certain small distance from the substrate are bound by 
the substrate. We estimate this distance as the pene- 
tration length C discussed in J2UJ. This is the length 
scale within which the polymers behave like unperturbed 
random walks rather than like strongly stretched chains. 
For a Na = 200 PMMA chains with radius of gyration 
R A w a (N A /6y/ 2 , we find ( » (R A /h A ) 4 ^h A « .5 h A . 
Then r can be found via 



1 



2D A 2 ho-hf 



(37) 




g{r)dr 




g{r)dr w 85%. (38) 
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In fact, this is an upper estimate for r, as density of 
ends g(r) has a singularity near the open end of the brush 
produced by the neglect of entropy of chain paths. The 
lower estimate for r would come from the assumption 
of constant g(r), so that r > ao/tiA ~ 10%, with ac- 
tual value being somewhere in between. Fortunately, this 
variation in r produces only an insignificant variation in 
the prefactor of D s ud- 

Finally, the self-diffusion constant for the chain flow 
along domain walls D c is also not readily available. 
Still, its value can be estimated from the analogy be- 
tween copolymer flow along domain walls and flow of 
star polymers, whose self-diffusion is given by D star w 
D arm exp(—N arm /N e ). Here N e « 180 is average num- 
ber of monomers per entanglement, and D arm is self- 
diffusion of linear polymers forming each arm. Then 
the self-diffusion coefficient for a PS star polymer with 
60,000 amu arms is Dps w .4 • 10~ 15 m 2 /s, and the 
self-diffusion coefficient for a PMMA star polymer with 
20,000 amu arms is £>pmma « A ■ 10" 15 m 2 /s*S Be- 
cause of similarity of these two coefficients, it is reason- 
able to assume that the self-diffusion coefficient for the 
whole block copolymer moving along domain walls is also 
D c w .4- 10~ 15 m 2 /s. Then the forced collective diffusion 
coefficient for the flow along bottommost terrace ()18fl is 
Dsiid ~ -7 • 10~ 15 m 2 /s, where known values for other 
parameters have been usect^S^i: ^ab ~ -6 • 10~ 3 N/m, 
h « 22 nm, and v ~ Na 3 « 1.7 • 10~ 25 m 3 . 

The final parameter influencing spreading dynamics 
(|2*5)l - (l3"6"j) is the activation free energy AE . Its mag- 
nitude strongly depends on the polymer conformation 
in the transition state. Two distinctly different confor- 
mations are possible 3 ^. In the first conformation, the 
hopping block propagates as a compact blob, thus mini- 
mizing its stretching energy at the expense of larger in- 
terfacial energy. In the second conformation, the hop- 
ping block is stretched between original and destina- 
tion domains, so that only a fraction of the hopping 
monomers comes in contact with the incompatible sta- 
tionary block. Using the calculations by Helfand 3 -, one 
can show that stretching mechanism dominates in the 
PS-PMMA films under discussion, and activation energy 
is roughly AE w xNpmmaJcT. 

Based on these quantities, we can now estimate the 
spreading rates of the PS-PMMA film. We made these 
estimates without any assumptions about the spreading 
pressure. This pressure depends on the molecular struc- 
ture and interfacial shape on an atomic scale at the ad- 
vancing front. It is difficult to estimate it without an ac- 
curate characterization of the surface and the monomers 
at this scale. Accordingly, we used the observed thick- 
ness at the front, namely hf instead of the pressure. It is 
observed to be about 0.85 times the equilibrium terrace 
thickness ho- Using this knowledge, we can apply Easl34l 
-E3 directly. 

First of all, the crossover time t c from the step-limited 
to the diffusion-limited dynamics, using Eq. 1371 is found 
to be about .25 hours. Hence diffusive flow along the 
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FIG. 6: Change of domain orientation during film thinning. 
The spreading of the film to the left carries a given region 
A in the step to the thinner region B, thereby reducing the 
angle 8 . 

terrace should dominate the spreading dynamics during 
most of the 24 hours of the experiment reported in Ref 0- 
Second, the typical distance L over which PS-PMMA film 
spreads in 24 hours at 523 K is predicted to be about 5 
microns (Eq. I34JI . This is indeed close to the observations 
by Hahm and Sibener-, where rim width was observed to 
increase by 1 - 3 microns. Finally, the spreading velocity 
after 24 hours is about 100 nm/hr. 

Domain alignment. It was already mentioned that 
hopping of chains between domains at terrace steps may 
account for alignment of the cylindrical domains in a 
spreading film. Specifically, in thin films, domain ori- 
entation is influenced by the nearby substrate, with do- 
mains usually oriented parallel or perpendicular to the 
substrate. For example, in PS-PMMA films on SiC>2 
lying-down cylinders are formed on nH— thick terraces 
(where H is domain period), while vertical cylinders are 
formed on (n + 1/2)7?— thick terraces. Spreading of the 
film and lateral growth of the terraces should then be ac- 
companied by changes in the observed domain patterns. 
Specifically, here we will discuss the growth and align- 
ment of the lying-down cylindrical domains on the grow- 
ing terraces. 

We first look at domain growth at the steps with hori- 
zontal domains on both sides of the step (Fig. [3J . This is 
the case, for example, at the "double steps" in the afore- 
mentioned PS-PMMA systems. Here term "double step" 
refers to the steps between nH— and (n + \)H— thick 
terraces. 

Following section IIVI topologically different domains 
of adjacent terraces may be assumed disconnected from 
each other (Fig. and the flow of material between ter- 
races is governed by single chain hopping events. The 
direction of this flow is determined by the difference in 
terrace chemical potentials, and during spreading goes 
towards stretched lower terrace. Then, upon arrival to 
the growing terrace, a copolymer chain has a strong en- 
ergetic bias to merge into existing domains rather than 
nucleate new ones. Hence the spreading should lead to 
the extension of the existing domains ensuring growth of 
continuous, defect free domain patterns. 

The thickness gradient in the step region, and to a 
smaller extent in the bulk of the terrace, can also account 



for the alignment of the growing surface domains in the 
thickness gradient direction. To demonstrate this align- 
ment effect, we first imagine that the cylindrical surface 
domains meet the transitional region at an oblique angle 
do to the growth direction, as shown in Figure [3] Since 
surface domains are not linked to the bulk domains of the 
thicker terrace, they have to end within the transitional 
region where the thickness is somewhat greater than that 
of the lower terrace. Then a small increment of motion 
of this system consists of two separate processes: first, 
a small amount of new material is added to the ends of 
the surface domains; second, these domains, being part 
of the growing lower terrace, move in the direction of the 
thickness gradient. Film thinning combined with conser- 
vation of volume then leads to a change in the direction 
of these domains, as described below. 

The growth increment extends the domains. This 
growth, we presume, occurs without domain bending so 
that #o remains unchanged, leaving the domain in the 
lowest energy unbent state. The translation of the step 
causes a given element of the film to compress in thick- 
ness and to lengthen in the direction of motion in order to 
conserve its volume (Fig. [SJ|. The resulting deformation 
resembles that of a piece of dough being rolled under a 
rolling pin. As domains thin, their angle Oq decreases. To 
quantify this decrease we notice that change in thickness 
of the film from ho to hf leads to its lateral extension 
Axj/Axq — ho/hf along the thickness gradient direc- 
tion x (Fig. EJ) , while lateral extension Ay along the step 
edge remains unchanged. Then the flattened domains 
emerge at an angle Of with the step edge, which is given 
by tan Of = Ayf/Axf = tanOo-ho/hf, and the extruded 
domains are bent towards the normal of the interface. 

This mechanism alone causes domains growing in an 
arbitrary direction to emerge partially aligned to the nor- 
mal. Further effects should cause this alignment to in- 
crease progressively with time. To demonstrate such an 
effect, we imagine that domains grow in increments of 
length s. During each of these increments domains ex- 
tend by distance s, and then propagate the same dis- 
tance along the gradient direction. The change in thick- 
ness that accompanies this motion, results in bending 
of domain tip from angle Oq to S with latter given by 

tan# s = tan#o • h s /ho = tan#o ^1 — ^^-s^j. The re- 
laxation of these bent domains into straight ones would 
require collective motion of large number of chains and 
exponentially long times. As this is highly unlikely, the 
next step of the growth will start with domain ends ori- 
ented at a smaller angle S to the normal, leading to 
progressive decrease in S with each growth step: 

Atan0 _ |V/i| 



and, hence, 

tan0 = tan O e~^ s . (39) 
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FIG. 7: Schematic top view of a step between terraces with 
lying-down semicylindrical domains on the growing lower ter- 
race and standing cylinders on the shrinking thicker terrace. 
Only the topmost terraces of the film are shown. 

Thus the growing domains progressively align in the 
normal to the terrace step direction. The characteristic 
length for this effect is determined by the size of the step 
region h/\Vh\ ~ 2 — 4 h so that characteristic length for 
alignment is several domain widths. 

To complete this section, we turn to an alternative 
combination of domain morphologies at terrace steps, 
with lying-down semicylindrical domains on the grow- 
ing lower terrace and standing cylinders on the shrinking 
thicker terrace (Fig. 0). In particular, this setup is of- 
ten observed in PS-PMMA films^ . While the general 
dynamics of the spreading should remain the same, some 
aspects of it may differ. First of all, the flow of material 
between terraces may proceed faster as some standing 
cylinders of the thicker terrace may continuously deform 
into short lying-down cylindrical domains (Fig. [7J that 
later merge with the existing domains of the lower ter- 
race. While this effect should have little effect on the net 
rate of spreading, which is primarily determined by col- 
lective diffusion in the bottommost terrace (not shown 
in Fig. UJ, it may lead to faster disappearance of the 
terraces with standing cylindrical domains, in agreement 
with experimental observations^. 

Another difference comes from the hexagonal pack- 
ing of the vertical domains, as it creates preferred di- 
rections for easy merging of vertical domains into hor- 
izontal domains. Since these directions do not depend 
on the wetting forces, there may emerge a misalignment 
between direction of spreading and direction of terrace 
growth, which was indeed observed in experiments on 
annuli structures^. 



VII. DISCUSSION 

We have shown that layered structure of spreading 
copolymer films should lead to different type of dynamics 
compared to macroscopic homopolymer or simple liquid 
drops. First of all, the spreading of copolymer films is 
determined by combination of two factors: the rate of 
chains flow inside constant-thickness terraces of the film 
and rate of chain flow between terraces. The hydrody- 
namic flow of chains along the terraces leads to the dif- 
fusion type law for spreading and its contribution 
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dominates at late times of spreading, with velocity de- 
creasing as v oc £ -1 / 2 . The hopping mechanism of chain 
flow at terrace steps, on the other hand, controls the 
early stages of spreading. The crossover time between 
the two regimes is given by (|37|) . It is notable that both 
of these regimes are different from the spreading observed 
in macroscopic drops, where spreading rate decreases as 
v oc t~ 9 . However, the v oc i -1 / 2 behavior is similar to 
the generic spreading rate of terraced fluids predicted be 
de Gennes and Cazabali 

Both regimes of spreading are predicted to exhibit 
strong dependence on chain polymerization. At early 
times (eq. I36fl , the dominant contribution to the spread- 
ing rate comes from the activation free energy AEq oc 
N a , so that v oc j\[- 2 - 7 e - N / N « , as opposed to v oc r\~ x oc 
A~ 3 in macroscopic homopolymer drops. Here a. = 1 for 
chains that hop between domains in a single blob, and 
a = 2/3 for chains that stretch between domains during 
hopping 33 . 

Dependence on N becomes even more complicated at 
large times, when diffusion along terraces dominates. It 
was mentioned above that the spreading rate at large 
times is primarily determined by chain flow in the bot- 
tommost terrace with collective diffusion constant D s nd 
(|18|l . For most copolymers, the main contribution is 
expected from the first term in l|18|l . which describes 
collective diffusion of fast chains through the matrix of 
substrate-bound slow chains. Then D s ud oc {h^rkd)^ 1 oc 
jy-4/3 e -A72Ar._ However, in systems with large degree of 
polymerization or small decrease in mobility £ near sub- 
strate, the second term becomes dominant with spread- 
ing rate independent of N: D s ud oc (hokg)^ 1 oc N°. 
Similarly, spreading rate should be independent of ./V if 
domains in the bottommost terrace obstruct continuous 
flow of chains, so that two-liquid approximation is inap- 
plicable and first term in l|18|) disappears. 

The flow of material between domains creates a dy- 
namic process of domain growth at the terrace steps. The 
chains entering the growing terrace preferentially merge 
the existing domains in that terrace, leading to their con- 
tinuous growth. In the previous section we have argued 
that this growth has preferential directions leading to the 
dynamic alignment of the domain. Before proceeding 
with further discussion of this model, it seems useful to 
compare it with other potential mechanisms for domain 
alignment, such as static, stress-induced alignments. 

In the annuli structures described above, stresses may 
come from the viscous flow, from the terrace steps, or 
may be caused by the change in circumference as material 
flows towards the center. In order to alter domain mor- 
phology, these stresses must be comparable to the inter- 
nal stresses in the film created by domain walls. In films 
with thickness h and surface tension between block jab , 
the internal stresses are roughly a wa ii ~ jAs/h ~ .6 • 10 5 
N/m 2 . 

This is significantly larger than the viscous stress 
produced by the interpenetrating brushes a V i SCOUS ~ 
r]*v/2<; ~ 2 • 10~ 2 N/m 2 (compare with l(T5)ll. making 



viscous stress an unlikely candidate for the cause of ob- 
served alignment. 

The other two stresses, on the other hand, come 
within one order of internal stresses, and could play 
some role in the alignment. Specifically, stress from a 
contact line is cr W etting ~ 7(^0 — hf)/ho ~ l^ 4 N/m 2 
(where hh and hf are thicknesses of relaxed and pulled 
films), and stress from the change in circumference is 
Oc ~ (j/h )(AR/R) ~ 10 4 N/m 2 (where R and AR 
are radius and width of the rim) . Despite the similarity 
of these stresses to the internal stresses in the film, the 
alignment caused by them would progressively improve 
towards the center of the rim, which contradicts the ob- 
servations on annuli structures^. 

The dynamic model of domain growth presented above, 
on the other hand, predicts better alignment on the 
grown parts of the terraces, that is in the central parts 
of the terraces and near the thicker terraces, in agree- 
ment with the experiment. This suggests that produced 
alignment is indeed caused by the dynamic effects accom- 
panying spreading. 

Our dynamical model of domain growth may provide a 
new way of producing aligned domain patterns. In partic- 
ular, any motion of terrace steps in the phase-separated 
film should lead to the growth of aligned domains on 
the growing terraces. One such experiment, besides the 
aforementioned experiment with annuli structures, is re- 
ported in Ref Il4l In this experiment, copolymer is de- 
posited onto a corrugated grating rather than a flat sub- 
strate. After the solvent has evaporated, the copolymer 
starts climbing up along the groove walls of the grating 
leaving aligned domain patterns inside the grooves. As 
this growth proceeds on time scales much longer than the 
thickness quantization, the observed alignment may be a 
result of the terrace growth. 

An alternative way to stimulate terrace growth may 
come via control over initial nucleation of islands and 
holes in the film. Specifically, under normal conditions, 
this nucleation happens at random places. However, if 
one applies periodic external catalyst of the island forma- 
tion, such as periodic temperature gradient, periodic elec- 
tric field, or, possibly, a standing sound wave, then the 
formation and growth of islands and holes may proceed 
in controlled manner leading to the growth of aligned 
domains. 



VIII. CONCLUSION 

We have shown that spreading of the phase separated 
copolymer films is different from the spreading of ho- 
mopolymer and other simple liquids. This difference 
is driven by domain walls inside copolymer and forma- 
tion of quantized thickness terraces. As a result, regular 
isotropic reptation mechanism is replaced by anisotropic 
flow of chains along domains, and slow activated flow be- 
tween domains. We have shown that at early stages, ac- 
tivated flow between domains at the terrace steps plays a 
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major role, giving constant velocity of spreading. At the 
late times, hydrodynamic flow inside terraces dominates, 
resulting in diffusion type velocity decrease v oc i -1 / 2 , 
similar to the generic spreading of terraced fluids pre- 
dicted be DeGennes and Cazabali 

The hopping nature and slow rate of chain flow at the 
terrace steps were shown to produce aligned and con- 
tinuous domain patterns on the growing lower terraces. 
Specifically, as new chains enter the growing terrace, they 
tend to merge with existing domains thus leading to their 
continuous extension. Furthermore, thickness gradient 
in the area of domain growth, i.e. near terrace steps, 
was shown to align growing domains perpendicular to 
the steps between terraces, thus leading to formation of 
aligned domain patterns on the growing terrace. 

This work shows that wetting flows offer a powerful 
means of controlling domains in copolymer films. The 
effect seems to arise not from flow-induced stress, but 



rather by flow-induced forced hopping of chains at ter- 
race boundaries. The alignment arises at these bound- 
aries, as in the conventional zone refining of crystals. As 
suggested above, one may imagine many ways to induce 
wetting flows in order to achieve a wide variety of mor- 
phologies. The potential for this new structuring mech- 
anism remains to be explored. 
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